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ABSTRACT 

We present a multi-wavelength, UV-to-radio analysis for a sample of massive (M* ~ 10 10 M Q ) 
IRAC- and MIPS 24/im-detected Lyman Break Galaxies (LBGs) with spectroscopic redshifts z^3 in 
the GOODS-North field (Ltjv>1-8xL*_ 3 ). For LBGs without individual 24/im detections, we employ 
stacking techniques at 24/im, 1.1mm and 1.4GHz, to construct the average UV-to-radio spectral 
energy distribution and find it to be consistent with that of a Luminous Infrared Galaxy (LIRG) with 
L IR =4.5^2'3 X 10 11 L Q and a specific star formation rate (SSFR) of 4.3 Gyr -1 that corresponds to a 
mass doubling time ^230 Myrs. On the other hand, when considering the 24/im-detected LBGs we 
find among them galaxies with Lir >10 12 Lq, indicating that the space density of z ~3 UV-selected 
Ultra-luminous Infrared Galaxies (ULIRGs) is ^(1.5±0.5) x 10~ 5 Mpc~ 3 . We compare measurements 
of star formation rates (SFRs) from data at different wavelengths and find that there is tight correlation 
(Kendall's r > 99.7%) and excellent agreement between the values derived from dust-corrected UV, 
mid-IR, mm and radio data for the whole range of Lir up to Lir ~ 10 13 L Q . This range is greater 
than that for which the correlation is known to hold at z~2, possibly due to the lack of significant 
contribution from PAHs to the 24/im flux at z ^3. The fact that this agreement is observed for 
galaxies with Lir > 10 12 L Q suggests that star-formation in UV-selected ULIRGs, as well as the bulk 
of star- formation activity at this redshift, is not embedded in optically thick regions as seen in local 
ULIRGs and submillimeter-selected galaxies at z = 2. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high redshift — infrared: 
galaxies 



1. INTRODUCTION 

One of the most fundamental quantities needed for un- 
derstanding the nature and evolution of galaxies is the 
star formation rate (SFR). To get reliable and meaning- 
ful estimates of the SFR for galaxies at high redshift, 
one needs a well defined sample of objects, coupled with 
multi-wavelength data that can provide a thorough and 
comprehensive investigation. 

One of the most successful methods of detecting high- 
z star-forming galaxies is the Lyman-break technique, 
pioneered by Steidel et al. (1996,2003). This technique 
has revealed a wealth of Lyman Break Galaxies (LBGs) 
at z ~3, now comprising an impressive catalogue of 
thousands star-forming galaxies at this redshift. Multi- 
wavelength studies of LBGs have provided extensive in- 
formation on various physical properties of these ob- 
jects. In particular, measurements at near-infrared wave- 
lengths and at 3.6-8/xm from the Spitzer Space Telescope 
IRAC instrument indicate that their stellar masses are 
typically 10 9 -10 n M Q (e.g., Shapley et al. 2001, Pa- 
povich et al. 2001, Magdis et al. 2010). 

The dust content and the SFR of LBGs at z ~ 3 are 
still poorly constrained. For their siblings at lower red- 
shift, z -2, Reddy & Steidel (2004) and Reddy et al. 
(2006), using multi- wavelength data ranging from X-rays 
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to radio, have reported that UV can be a reliable SFR 
indicator if corrected for dust attenuation by an aver- 
age factor between 4.4 and 5.1. The validity of the UV 
as a robust SFR indicator has also been presented by 
Daddi et al. (2005, 2007) for a sample of near-infrared se- 
lected galaxies at 1.5 < z < 2.5 identified using the BzK 
technique (see also Dannerbauer et al. 2006). A similar 
multi-wavelength study for the z ~3 LBGs, though, is 
still needed. 

In this letter we make use of the unique compilation of 
multi-wavelength data on the Great Observatories Ori- 
gins Deep Survey North field (GOODS-N) to explore the 
SFR and the infrared luminosities (Lir) of z ~3 LBGs. 
Our aim is to fully characterize the spectral energy dis- 
tribution (SED) of a typical LBG from rest-frame UV 
to radio wavelengths, to compare different tracers of star 
formation, and to test whether the UV can provide a re- 
liable measurement of star formation at z ~3. For this 
letter we adopt a A cold dark matter (ACDM) cosmology 
with H = 71 km s^ 1 Mpc" 1 , Vt m = 0.27 and fl A = 0.73, 
while the magnitudes presented in this work are all in 
the AB magnitude system. 

2. THE LBG SAMPLE AND DATA SETS 

2.1. Rest-frame UV data 

In this study we selected 69 LBGs in the northern 
field of the Great Observatories Origins Deep Survey 
(GOODS). Their original selection was based on their 
optical colours (U n ,G,R) to R =25.5 by Steidel et al. 
(2003) and subsequent optical spectroscopy (Steidel et 
al. 2003, Reddy et al. 2005) has confirmed their high 
redshift nature with a median z=2.95. Optical spec- 
troscopy has also been used to determine the absence 
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Fig. 1. — Stacked image at 24/im (a), 1.1mm (b) and radio 1.4GHz (c) of IRAC-LBGs d) Stacking simulations at 24/xm. Distribution 
of the measured fluxes derived from 50.000 stackings at 40 random positions along with the best gaussian fit (rms=1.04/ijy). The red line 
denotes the flux measurement of the stacking at the position IRAC-LBGs. This indicates a ~7.9<r detection at the stacked position of the 
LBGs. e) Stacking at 1.1mm. Same as in d) but for the Aztec 1.1mm and for 5000 stackings. This figure indicates a ~3.7<r detection at 
the stacked position of the IRAC-LBGs. 



of AGN signatures (i.e., strong high ionization emission 
lines) in their rest-frame UV spectrum, although a deeply 
obscured AGN cannot be ruled out from these data. We 
also use BViz data obtained from the Advanced Camera 
for Surveys (ACS) on board the Hubble Space Telescope 
(HST) (Giavalisco et al. 2004). 

2.2. IRAC and MIPS 2^m data 

Rest-frame NIR identification and IRAC photometry 
(5(7, fs. 6=0. 21 /iJy) of the LBGs in our sample have been 
presented by Magdis et al. (2008,2010), who also showed 
that their mid-IR colours are consistent with those of 
star-forming galaxies at z ^3. Here, we focus on the 49 
LBGs with at least one IRAC detection ([3.6]ab <25.0), 
in order to place a lower mass limit of M* ^1O 1O M 
(Magdis et al. 2010) and facilitate a robust investiga- 
tion of their properties and their average SED. We not 
the R magnitude limit of our sample, (Rab <25.5), cor- 
responds to Lyv > l-8xL* =3 (Mj^=-21, Steidel et al. 
1999) 

We matched our sample with a GOODS-N 24/im cat- 
alog (5cr ~ 20^Jy, translated to SFR«350 I^yr" 1 or 
Lib, «2x10 12 L q based on Chary & Elbaz 2001 mod- 
els at z=3) produced by the GOODS team (Dickinson 
et al. in prep.), and searched for counterparts within a 
2" diameter separation centered on the optical position. 
We identify 9 LBGs and we add another five 24/im- 
detected LBGs in the Extended Growth Strip (EGS, 
5cr ^70 /xJy) by Rigopoulou et al. (2006), to increase 
our 24/im-detected sample. We note that the original 
optical selection of these extra five EGS MlPS-selected 
LBGs is identical to that of the GOODS-N sample, so 



they share similar UV properties and including them in 
our sample doesn't introduce any bias. Henceforth, we 
will refer to LBGs that are individually detected at 24/im 
as MIPS-LBGs (14 objects), and those that are not as 
IRAC-LBGs (40 objects). In practice, all MIPS-LBGs 
are also detected in all four IRAC bands. The redshift 
range of the two sample is 2.63< z < 3.41 (IRAC-LBGs) 
and 2.60< z < 3.31 (MIPS-LBGs) following a similar dis- 
tribution. The median redshift of the MIPS- and IRAC- 
LBGs is 2.92 and 2.98 respectively. 

For the 24/mi-undetected IRAC-LBGs in GOODS-N 
we employed median stacking analysis. We first sub- 
tracted all detected sources using the PSF used for the 
source extraction, and cut sub-images centered at the op- 
tical position of each undetected LBG. To avoid contam- 
inating the stacked signal from residuals, we only added 
galaxies to the stack if there were no bright MIPS sources 
within ^4" of those galaxies. Then a stacked flux was 
measured in a manner similar to the measurement of 
the detected MIPS sources. The final stacked image at 
the position of the IRAC-LBGs is shown in Figure la. 
To quantify the error of our measurement we stacked at 
40 random positions and repeated it 50.000 times. As 
expected for white noise, the distribution of the fluxes 
follows a gaussian shape with an rms of 1.04 fiJy (Figure 
Id) that we adopt as the uncertainty of our measure- 
ment. The median flux density of the IRAC-LBGs as 
derived from stacking is / 2 4=7.91±1.04/zJy (S/N -8). 

Previous studies have demonstrated that MIPS-LBGs 
are, on average, more massive, relatively older and 
dustier compared to LBGs that are undetected at 24/im 
(e.g., Rigopoulou et al. 2006, Magdis et al. 2010). In the 
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rest of this paper we study the IRAC-LBGs and MIPS- 
LBGs separately to enable a comparison between the two 
sub-populations of LBGs and examine our results as a 
function of Lir and dust extinction. For our statistical 
analysis we always refer to median values. 

2.3. Aztec 1.1mm data 

Recently, a deep (a ~ 0.96-1.16 mJy beam -1 ) and uni- 
form 1.1mm survey of the GOODS-N field with AzTEC 
(Wilson et al. 2008) was conducted on the James Clerk 
Maxwell Telescope. Matching our sample with the cata- 
logs published by Perera et al. (2008) and Chapin et al. 
(2009), returned no individual detection down to 3.75ct 
level. Since none of our LBGs is individually detected 
we used the publicly available maps to do stacking. We 
perform median stacking at the positions of IRAC-LBGs 
after rejecting two of them that are separated by less 
than an Aztec beam half-width (9") from any known 
Aztec detection (Figure lb). The median flux that we 
recover for the IRAC-LBGs is /i.i mm =0.41±0.11 mJy 
(3.7a). The robustness of this flux was also tested by 
the same method employed for the 24/zm stacking (Fig- 
ure le) and bootstrapping. Finally, stacking at the po- 
sitions of the seven GOODS-N MIPS-LBGs (two were 
rejected due confusion) returned no detection indicating 
a 4cr upper flux density limit of 1.08 mJy. 

2.4. VLA 14GHz radio data 

We use a new deep 20 cm (1.4 GHz) imaging of the 
GOODS-N field obtained at the VLA in ABCD config- 
urations (G. Morrison et al. 2010, submitted). Previ- 
ous observations by Richards (2000) contained 40 usable 
hours of A-array data which was supplemented with 125- 
hr yielding an rms of 3.9^Jy beam -1 near the phase 
center. We matched our LBG sample to the 5a radio 
catalog (Morrison et al. 2010) and identified one MIPS- 
LBG, HDFN-M23 with a flux Si 4 gh z =21.19±4.2 (j,Jy 
(2=3.21). 

In order to reach deeper radio flux densities, we stacked 
individually undetected LBG sources using the following 
technique. Using a primary beam corrected radio im- 
age, we made lOOpixel x lOOpixel sub-images centered 
on the LBG positions. The exported FITS files were then 
stacked using the Terapix co-add software Swarp (Bcrtin 
et al. 2002). The stacking used the pixel reference frame 
and the resulting stacked image was median combined. 
We stacked separately MIPS- and IRAC-LBGs and de- 
rived a median flux density of 8. 5 ±2. 2 /iJy and 3. 6 ±0.8 
^Jy respectively. The stacked image of the IRAC-LBGs 
is shown in Figure lc. 

3. SFR INDICATORS 

UV SFR estimates. To determine the UV-corrected 
SFR (SFRuvcor), we use the GOODS ACS photometric 
catalog retrieved from MAST STScI (Version v2.0). At 
z ^3 the observed V and I and z bands correspond to 
rest-frame 1500 A, 2000 A and 2400 A respectively, al- 
lowing a robust estimate of the f3 slope. We use dust free 
models of continuous SFR, solar metallicity, age t s f r = 
lOOMyrs and Salpeter IMF, generated with the new code 
of Chariot & Bruzual 2007 (CB07 private com.), to fit 
the SED of each individual LBG after correcting for red- 
dening using the Calzetti (2000) attenuation law and 



correcting for the IGM attenuation using the prescrip- 
tion of Madau (1995). From the best fit model, we de- 
rive E(B-V) values and apply a K-correction to infer the 
observed and subsequently the intrinsic flux density at 
1500A in the rest frame. The average E(B-V) value is 
0.16 (Ay=4.6) and agrees well with previous studies (i.e., 
Shapley et al. 2001). The adopted CB07 models yield a 
relation between the SFR and the monochromatic 1500A 
luminosity given by: 

SFR(M yr- 1 ) = L^erg s^Hz" 1 ]/^ x 10 27 ) 

(1) 

Using extinction-corrected L ^ in equation 1, provides 

an estimate of the total SFR (SFRuv, corrected)- The to- 
tal SFR can be considered as the sum of an unobscured 
component (SFRuv, uncorrected), which can be computed 
from equation 1 by using the observed L °, and an 

1 1500A' 

obscured component (SFR ODS c), corresponding to the en- 
ergy absorbed by dust. This absorbed energy is in turn 
rcradiated at mid- and far-infrared wavelengths, SFR b sc 
= SFRir, and it is this latter quantity that we will infer 
from the mid- and far-infrared data. Therefore, we may 
write: 

SFRtjV, corrected = SFRtr + SFRuv, uncorrected (2) 

24/mi, 1.1mm, 1.4GHz radio SFR and Lir esti- 
mates. To convert 24/zm and Aztec fluxes to total (8- 
1000/mi) Lir, we use the luminosity-dependent SED li- 
brary of Chary & Elbaz (2001) (CE01), while for com- 
parison we also consider the Dale & Hclou 2002 (DH02) 
models (with luminosities normalized as described by 
Marcillac et al. 2006) and SED templates of Arp220 and 
M82. This is done by interpolating Lir over the tem- 
plate SEDs, sorting the value of Lir that corresponds to 
the observed 24 /im flux density. For the radio-based Lir 
estimates, we first derive the observed radio luminosities 
and then assume a radio spectral index of a=-0.8 to get 
the rest-frame 1.4GHz radio luminosities. Then we use 
the local radio-IR correlation (Condon 1992) to deter- 
mine Lir : 

L m /L Q = 3.5 x 10~ 12 L(1.4 GHz) [WHz^ 1 ] (3) 

The Lir is subsequently converted to SFRir using 
Kennicutt (1998) : 

SFR radio [M yr- 1 ] = 1.73 x lO- lo L m [L ]. (4) 

4. COMPARISON OF UV, MID-IR, MM AND RADIO 
SFR/L m ESTIMATES 

4.1. UV vs 24/j.m 

We first compare the SFR and Lir derived from UV 
and mid-IR in Figure 2, where we also show the expected 
correlation based on the CE01, DH02, Arp220 and M82 
templates. For the IRAC-LBGs, we use the Lir derived 
from the stacking analysis, while for the rest we use the 
one derived from their measured 24/im flux densities. It 
is evident that between UV- and 24/im derived luminosi- 
ties based on CE01 and DH02 templates, there is a close 
agreement which for the case of the latter templates is 
more prominent. The Kendall's r test detects a corre- 
lation at a >99.7% confidence level, that holds up to 
Lir ^10 13 L Q , where the two estimates appear to de- 
viate. The two galaxies with the largest 24/xm-derived 
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Fig. 2. — A comparison between the 24/im inferred luminosi- 
ties and those derived from the obscured UV-light. The blue cir- 
cle shows the stacking result for IRAC-LBGs, black squares show 
LBGs with 24/im detection (MIPS-LBGs) while red squares show 
MIPS-LBGs from EGS which have a 1.2mm MAMBO detection. 
Blue arrows indicate 24/im upper limits (3rr) for the IRAC-LBGs. 
For these symbols, y-axis indicates Ljr (left) and SFR, both de- 
rived from UV uncorrected measurements. The black line shows 
the one-to-one relation based on CE01. The colored lines show the 
correlation between CE01 and DH02, ARP220 and M82. For these 
lines the left vertical axis should be read as the Ltr predicted from 
the corresponding models. 

Lir and with the largest deviations between the UV 
and 24/im-derived SFR, are also detected at 1.2mm by 
MAMBO (Rigopoulou et al., in prep.), and are therefore 
members of the SMG population for which such a trend 
is known to exist either due to star formation embedded 
in optically thick regions or due to contribution from an 
AGN to the mid-IR output, or both, (e.g., Chapman et 
al. 2005, 2009, Pope et al. 2008). We also note that for 
these two galaxies, a scaled-up M82 template provides a 
good agreement between the UV- and 24/im derived Lir. 

A similar correlation between UV and mid-IR derived 
SFRs has also been observed at z ^2 for BzK (Daddi et 
al. 2007) and BX/BM (Reddy et al. 2006) galaxies with 
the two estimates though deviating at lower luminosities 
(^3xl0 12 L©) when compared to our findings at z^3. 
It is therefore indicated that the UV-24/im correlation 
at z~3 is better than that found at z^2. One possible 
explanation could reside in the contribution of PAH fea- 
tures in the observed mid-IR flux. While evidence was 
found that PAHs are enhanced in z^2 star forming galax- 
ies (Murphy et al. 2009) with respect to local galaxies, 
this may not be the case anymore at z~3 when galaxies 
were less metal rich and PAHs less abundant. We note 
that while for z>~2 galaxies the 24/im is centered at the 
PAH emission, at ~3 the band traces only the 6.2/im and 
part of the 7.7/im PAH feature. This shift though, can- 
not fully explain our findings as for a M82-like object, the 
contribution of the PAHs to the the 24/im flux density is 
comparable at the two redshifts (50-70%). Alternatively, 
it could be argued that the mid-IR radiation of our sam- 
ple is less polluted by the emission of hot dust heated 
by an AGN, which could be a natural result of the rapid 



drop of the number density of AGNs above z=2 (Wall et 
al. 2005). 

4.2. UV vs radio and 1.1mm 

One independent way to check the UV SFR estimates 
is to compare them against radio observations. In Figure 
3, we plot the SFR estimates based on the radio fluxes 
and the UV corrected for dust extinction for three sam- 
ples: HDFN-M23 which is individually detected in the 
radio map, the stacked flux of radio-undetected MIPS- 
LBGs, and the stacked flux of radio-undetected IRAC- 
LBGs. We see that there is an excellent agreement be- 
tween the radio and UV estimates, testifying to the va- 
lidity of UV as a SFR indicator for UV-selected galaxies 
at z ~3. The corresponding SFRs derived from UV (and 
radio) are 90t^ M yr- 1 ,(96±32 Moyr" 1 ) for the IRAC- 
LBGs, 250±|^ M yr-\(28O±85 M^r" 1 ) for the MIPS- 
LBGS and 808 M Q yr-\(870±200 M^r" 1 ) for HDFN- 
M23. A similar study by Carilli et al. (2008), comparing 
radio and UV data for z ~ 3 LBGs in COSMOS, inferred 
an average UV attenuation factor of ~ 1.8, smaller than 
the average value of ~ 5 derived here and in other stud- 
ies (e.g., Reddy & Steidel 2004; Reddy et al. 2006). We 
argue that the discrepancy arises from the fact that we 
focus on massive (dusty) LBGs with robust spectroscopic 
redshifts and compute dust attenuation on an object-by- 
object basis from the UV spectral slopes, whereas Carilli 
et al. primarily use photometric redshifts, infer only an 
average attenuation and do not consider a mass-limited 
sample. Finally, we explore the relation between the UV 
and the Lir as derived from mm. Following the same 
prescription as above, we convert the mm stacked flux 
the IRAC-LBGs and find it to be in agreement with the 
Lir derived from the UV. The two corresponding val- 
ues are (6.2±3)xlO n L Q (mm) and 5.1 xlO 11 L Q (UV). 
Converting the stacked mm value to SCUBA850/im by 
using the formula presented by Ivison et al. (2005) we get 
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Fig. 3. — A comparison between the radio 1.4GHz inferred 
SFR and those derived from the UV light from star-formation 
corrected for obscuration. Blue square shows the radio stacking of 
the IRAC-LBGs, black square shows radio stacking of MIPS-LBGs 
while green square shows estimates for HDFN-M23, an individu- 
ally detected LBG in the radio map. Straight line shows one to 
one correlation. The arrows indicate upper limits for the IRAC- 
(blue) and MlPS-detected LBGs. 
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Fig. 4. — Average SED of a typical IRAC-detected LBG at z ~3. For the SED we use the median BViJK+IRAC photometry of 
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S , 850 =0.85±0.27mJy indicating that IRAC-LBGs emit at 
the sub-mJy level at 850^m. We note that similarly 
to the SMGs, where UV is an unreliable indicator of 
SFR, there is evidence of existing optically thick star- 
forming regions in the case of sub-mm detected LBGs 
(e.g., Rigopoulou et al. 2010 in prep, Chapman & Casie 
2009). 

5. DISCUSSION 

Considering the median fluxes of the IRAC-LBGs for 
the rest-frame UV to NIR (i.e U„ to 8.0/im) and stacked 
fluxes at 24/im, l.lmm and 1.4GHz, we construct the av- 
erage SED of a typical (24^m faint) IRAC-LBG. We fit 
the rest-frame UV to NIR with model SEDs generated 
using the CB07 code, and the mid-IR to radio with CE01 
templates. The photometric points along with the best 
fit model are shown in Figure 4. The best-fit CB07 model 
indicates an average stellar mass of M* ^2.2xl0 10 M Q 
and an average SFR ^85 M Q yr _1 . Using this mass esti- 
mate and the SFR derived from our multi-wavelength 
analysis we derive a specific SFR (SSFR, defined as 
SFR/M Q ) ~4.3 Gyr" 1 , corresponding to a mass dou- 
bling time of ~ 230Myrs. This value is very close to the 
one presented by Magdis et al. (2010) (4.5 Gyr^ 1 ) and is 
larger than that found at lower and higher redshifts, re- 
inforcing their argument that the evolution of the SSFR 
peaks at z ~3. Based on CE01 models we find that the 
average IRAC-LBG is a LIRG with Li R =4.5±^xl0 11 
Lq. We note that Lir corresponds to a typical dust 



temperature Td=35K, based on Td measurements of lo- 
cal LIRGs (Yang et al. 2007). 

On the other hand, MIPS-LBGs have higher lumi- 
nosities, indicating that ULIRGs are present among the 
UV-selected galaxies at z ~3. For instance, HDFN- 
M23 which was individually detected in the radio, has 
Lir ~5(±2)xl0 12 L©. We calculate the comoving vol- 
ume for the redshift range of our sample (2.5< z <3.5) 
and find that the space density of ULIRGs (based on UV 
Lir) LBGs is ~1.5±0.5xl0~' 5 Mpc~ 3 . This is a factor 
of ~10 smaller than the space density of the z ~2 BzK- 
selected ULIRGs (Daddi et al. 2007). We note that the 
UV selection is likely to miss more obscured, UV faint 
ULIRGs as our sample is limited to Rab < 25.5 while, 
for example, an Arp220 like object at z ~3 would have 
Rab -26.9. 

The good agreement between SFR indicators that are 
affected by dust extinction (UV) or not (MIR, mm, ra- 
dio) suggests that the bulk of the star formation activity 
in massive UV-selected galaxies takes place in optically 
thin regions. Since by focusing on the massive LBGs we 
also select those that are most affected by dust extinction 
(e.g., Magdis et al. 2010), it is reasonable to extend this 
result to the whole population of UV galaxies including 
the less massive, less dusty LBGs. Furthermore, the fact 
that the agreement holds even for the case of the most 
massive, dusty, 24/im-detected LBGs with Lir >10 12 Lq, 
indicates that UV-selected ULIRGs at z=3 are transpar- 
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ent to UV light, meaning that we can estimate its SFR 
given its rest-frame UV 1500 A luminosity and UV slope, 
contrary to the local ULIRGs and z=2 SMGs. Similar 
results have been reached by Reddy et al. (2006) for 
UV-selected and by Daddi et al. (2007) for BzK selected 
z~2 LIRGs and ULIRGs. Both studies find that galaxies 
of a given bolometric luminosity are on average a factor 
of 8-10 less dust obscured at z^2 than at the present 
epoch. The fact that our study suggests that at z^3 
ULIRGs are optically thin at even higher Lir than that 
at z^2, indicates even less obscuration at z~3 confirm- 
ing the trend between galaxies at z = 0, z ~ 1, and z 
~ 3 (Adelberger & Steidel 2000). This could plausibly 
be explained as a result of increasing dust-to-gas ratios 
as we move from the high-z to the present universe. As 
galaxies evolve, they convert gas into stars which in turn 
enrich the interstellar medium with dust. If the dust 
distribution becomes more compact with time (assum- 
ing mergers that drive dust and gas to the center of the 
galaxy) the overall result would be an increase of the 
dust column density toward star-forming regions making 
ULIRGs progressively optically opaque at later epochs. 

To put constraints on the census of the cosmic star 
formation at z=3, we should also consider the missing 
fraction of star formation embedded in optically thick 



regions, that takes place in galaxies not selected in UV. 
On the other hand, there is evidence of decreasing ob- 
scuration with increasing redshift for a given Lir (e.g., 
Reddy et 2008), pointing towards less optically thick 
star-formation at higher redshifts. Combining the above 
with the fact that the bulk of the currently known SMGs 
are at z ^2.2 (Chapman et al. 2005), we can assume that 
at z=3 the contribution of the SMGs to the total SFR 
density, is not dominant, and hence locate the census of 
high-z star-formation in optically thin regions. We stress, 
that such a scenario cannot be confirmed based solely on 
our sample, as it has been shown that the Lyman Break 
technique can miss a large fraction of massive (dusty) 
galaxies at high-z (e.g., Daddi et al. 2004, Van Dokkum 
et al. 2006). 

Further insights into the far-IR properties of the LBGs, 
such as dust temperature and dust mass will be provided 
by deep surveys with the Photodetector Array Camera 
& Spectrometer (PACS, 70-,100-,160^m) and the Pho- 
tometric Imaging REceiver (SPIRE, 250-,350-500^m) on 
board the Herschel Space Observatory (HSO). Based on 
the average spectrum of Figure 4, the predicted flux 
densities in the HSO bands of a typical IRAC-detected 
LBG are /i O o=0.13 mJy, / 160 =0.47 mJy, / 25 o=l-18 mJy, 
/ 350 =1.52 mJy and / 500 =1.22 mJy. 
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